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We express the high-frequency (anti-adiabatic) hmit of the exchange-correlation kernels of an in- 
homogeneous electron gas in terms of the following equilibrium properties: the ground-state density, 
the kinetic stress tensor, the pair-correlation function, and the ground-state exchange-correlation 
potential. Of these quantities, the first three are amenable to exact evaluation by Quantum Monte 
Carlo methods, while the last can be obtained from the inversion of the Kohn-Sham equation for 
the ground-state orbitals. The exact scalar kernel, in this limit, is found to be of very long range 
in space, at variance with the kernel that is used in the standard local density approximation. The 
anti-adiabatic xc kernels should be useful in calculations of excitation energies by time-dependent 
DFT in atoms, molecules, and solids, and provides a solid basis for interpolation between the low- 
and high-frequency limits of the xc kernels. 



I. INTRODUCTION 

The exchange-correlation (xc) kernel is a quantity of 
central importance both in time-dependent density func- 
tional theory (TDDFT) and in time-dependent current 
density functional theory (TDCDFT). It is formally de- 
fined as the functional derivative of the time-dependent 
xc potential Vxc with respect to time-dependent den- 
sity in TDDFT, or as the functional derivative of time- 
dependent xc vector potential Axc,i with respect to time- 
dependent current-density jj in TDCDFT. In practice, 
this kernel connects the physical density- and current re- 
sponse functions of the interacting many-body system to 
those of a fictitious non-interacting system - the so-called 
Kohn-Sham system - which has the exact ground-state 
density. More precisely, one has 
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in TDCDFTi^ Here x and x denote the exact density and 
current response functions (the second being a tensor, 
as emphasized by the hat), Xs a-nd Xs ^-re their Kohn- 
Sham counterparts, uj is the frequency, c is the velocity 
of light, e is the absolute value of the charge of electron, 
and Vc{r) = e^/r is the Coulomb interaction. 

The great importance of the fxc kernels stems from 
the fact that the Kohn-Sham response functions usually 



misrepresent both the response and the excitation spec- 
trum of the system. Poles of Xs, for example, occur at 
differences of Kohn-Sham eigenvalues, which are known 
to be, at least conceptually, unrelated to the true exci- 
tation energies of the system. It is the task of fxc to 
shift the poles of the KS response functions from these 
unphysical values to the actual physical ones, and this is 
a very difficult task, particularly when it comes to band 
gaps and excitons. 

Very few exact results are known about fxc- Even in 
the simplest case of a homogeneous electron gas (in which 
case fxc and fxc are functions of |r — r'| and fxc decom- 
poses into independent longitudinal and transversal com- 
ponents) very little is known. The static limit of fxcili ^) 
has been calculated by variational and diffusion Quantum 
Monte Carlo (QMC) and fitted to simple formulasj^ii It 
provides the basis for the so-called adiabatic approxima- 
tion, which is based on the assumption that the exchange- 
correlation potential relaxes to the ground-state form on 
a time scale that is much shorter than the time scale of 
the evolution of the system. In the opposite limit of high- 
frequency, the quantity lim„_j.oo fxd'-^) = f!^ is known 
from exact sum rules to be expressible in terms of the 
density, the ground-state kinetic energy, and the ground- 
state pair correlation function (or equivalently the static 
structure factor) , all of which can be obtained from QMC 
calculations. This limit provides the basis for an "anti- 
adiabatic approximation" , in which fxc is interpreted as 
the stiffness of the ground-state wave function with re- 
spect to instantaneous virtual deformations. This point 
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will be expounded more fully in the following. 

The frequency dependence of the homogeneous kernels 
has been the object of several paramctrizations, begin- 
ning with the one by Gross and Kohn in 1985^ (see also 
Refs. 0and@), which interpolate between low- and high- 
frequency limit, and incorporate as many exact proper- 
ties and constraints as possible. Perturbative calcula- 
tions, which treat electron exchange exactly, have also 
been carried out.— 

Historically, the xc kernels of the homogeneous elec- 
tron gas have played (and still play) an important role as 
the basis for local density approximations, whereby the 
xc kernels of the non-homogeneous systems are approxi- 
mated in terms of those of a homogeneous system, eval- 
uated at the local density. However, it has become clear 
in recent years that some features of the exact kernels 
cannot be reproduced by the local density approxima- 
tion. For example, it is now known that the exact scalar 
xc kernel has a singularity at finite frequency and q ^ 
(corresponding to long range in space), which plays an 
important role in the optical spectra of extended system. 
Reining et al^ introduced phenomenologically a singu- 
lar fxc, and using it "on top" of bands calculated by the 
GW method demonstrated considerable improvement in 
the calculated optical spectrum of silicon. One way to 
introduce a singular f^c from first principles is to derive 
it from a local density approximation for the tensor f^c 
of TDCDFT. Indeed, any regular approximation for f^c, 
when translated into an approximation for f^c, does pro- 
duce a singularity of the expected formJ^iii The extent 
to which this may help to resolve the difficulties in the 
calculation of optical spectra remains unclear. Further it 
now appears that not only the singular behavior of fxc at 
q = 0, but also the detailed non-local spatial dependence 
of this kernel plays a role in the calculation of the optical 
spectra. 

In view of these difficulties it seems particularly timely 
to try and learn as much as possible about the exact prop- 
erties of fxc in non-uniform systems. In a recent paper 
we have taken a first step in this direction by deriving 
an exact expression for the fxc of a weakly inhomoge- 
neous system to second-order in the strength of the ex- 
ternal potential This expression involves the kernels of 
the homogeneous electron gas (which, as discussed above, 
are still poorly described), but is otherwise fully nonlocal. 
Unfortunately, the perturbative approach restricts us to 
metals, and seems to exclude a priori the interesting case 
of semiconductors. 

In this paper we calculate the exchange-correlation ker- 
nels for arbitrary non-homogeneous systems in the anti- 
adiabatic limit. As remarked above, the anti-adiabatic 
limit of the homogeneous exchange-correlation kernel is 
well known (a state-of-the-art formula for it is provided 
in the appendix). Here we extend that knowledge to the 
considerably more complex case of a non-homogeneous 
system. In spite of the additional complication, we are 
still able to express the result in terms of ground-state 
density, kinetic stress tensor, pair correlation function 



and exact exchange-correlation potential. All these quan- 
tities are well within the reach of modern numerical al- 
gorithms, and therefore our results open the way to an 
accurate study of fxc in the anti-adiabatic limit. 

The high-frequency response of quantum many-body 
systems has recently come under intense scrutiny in the 
context of the continuum mechanics formulation of quan- 
tum many-body theory.— In this theory one attempts to 
describe the dynamics in terms of a displacement field 
u(r, i), describing the displacement of an infinitesimal 
volume element from its "equilibrium position" in the 
ground-state (in this approach the time-evolution is de- 
scribed as a geometric deformation of the ground-state). 
It turns out that the simplest approximation that leads 
to an eigenvalue problem for the excitation energies is 
in fact equivalent to a high-frequency expansion for the 
current response function. We get 

mu(r,t) = -VVi{r,t) 

^— [ dr'M(r,r')-u(r',t), 

"-o(r) J 

where no(r) is the ground-state density, and the tensor 
M(r, r') is the first moment of the current excitation 
spectrum [see Eq. ([5])]. This quantity can be expressed 
in terms of ground-state properties, i.e., it is the second 
variational derivative of the ground-state energy with re- 
spect to the instantaneous displacement field u(r), eval- 
uated at u = 0. We will see that M(r, r') also plays a 
crucial role in the construction of . 

This paper is organized as follows. In Section |ll] we 
present a general formulation for the anti-adiabatic ten- 
sor and scalar xc kernels in terms of the first spectral 
moment of the current density fluctuation spectrum. In 
Section Hill we provide explicit expressions for these ker- 
nels in finite or infinite systems in real space. In Sec- 
tion IIVI we make connection with the physically instruc- 
tive continuum mechanics formulation, and discuss the 
interpretation of fxc in terms of elastic constant (stiff- 
nesses) of the electronic system. In Section |V] we present 
the expressions for the tensor and scalar xc kernels for 
an infinite periodic system in momentum space. Work- 
ing in momentum space allows us to focus most clearly 
on the long-rangcdness of the kernels in the g — ?> limit: 
Section [VII summarizes our results for the singular parts 
of the kernel. Finally, Appendix |X] presents a detailed 
derivation of the all- important double commutator M, 
and Appendix [B] presents results for the longitudinal and 
transversal component of the anti-adiabatic fxc in the ho- 
mogeneous electron gas, obtained by using state-of-the- 
art values of the static structure factor and correlation 
energy. 



3 



II. GENERAL FORMULATION 



The current response tensor can be written as 
, eno(r)(5(r-r') „ 



cm 



+e^[F(a)-n/3)] 

a,/3 



{a\j,{rm{(3\Mr')\a) 
LJ - Ep + Ea + iri 



, (3) 



where |a) and Ea are the many-body wave functions and 
eigencnergies, respectively, no(r) is the ground-state elec- 
tron density, j(r) = X;„[Pn^(r - r„) (5(r - r„)p„]/(2m) 
is the paramagnetic current-density operator, p„ is the 
momentum operator of the n-th electron, -^(q;) are the 
diagonal elements of the equilibrium (or ground-state at 
T = 0) density matrix, w is the frequency, c is the ve- 
locity of light, e and to are the absolute value of the 
charge and mass of electron, respectively. Expansion of 
the right-hand side of Eq. ([3]) in l/w gives to the second 
order 



X^i.(r,r',a;) 
where 



uj—^oo 6 
C 



,(4) 



(5) 



is the first moment of the imaginary part of Xpi'i the 
angle brackets denoting the ground-state or equilibrium 
average at zero or finite temperature, respectively. 



2m 



n '-ml 



is the Hamiltonian of the unperturbed system with Vo(r) 
being the external potential. Inversion of Eq. yields 



tj— >oo C 

e 



m6{r-r') A'V(r,r') 



«o(r) 



uj'^no{r)no{r') 



Using the definition of the tensor exchange-correlation 
kernel of TDCDFT? of Eq. Q we can write 



ea;^ri,o(r)rio(r') ' 



(6) 



where 



Af;,= (r,r') = M^,(r,r')-M;,(r,r') 

- noir)no{r')d^dlVc{\r^r'\). 



The connection to the scalar quantities follows from 
the equality 



(8) 



where x(r, r',w) is the density response function, and 
the summation over the repeated coordinate indices is 
implied. By virtue of Eqs. (jlj and ([8|) we can write 



X{r,r',u) — ^9^no(r)a|^J(r - r') 



the inversion of which gives 



(9) 



X (r,r ,w) « moj g{r,r) 
- I [d'^Gir, r")]AV(r", r"')[9;"e(r", r"')]dr"dr"', 

(10) 

where Q{r,r') = [dfj,no{r)d'^5{r — r')]~^ can be written as 
the solution of the following equation 

-d^noir)d^g{r,r') = div-r'). (11) 

With the use of the definition of the scalar xc kernel 
of Eq. (P we have from Eqs. ^ and 

/,°^(r, r') = J [d';g{r, r")]A./^-(r",r"') [^^(r", r"')]dr"dr"' , 

(12) 

where /^(r,r') = lim^^oo ^^(r, r', w). 

III. TENSOR AND SCALAR XC KERNELS IN 
REAL SPACE 

In Appendix]^ we outline a rather lengthy evaluation 
of the double commutator in Eq. ([5]) which leads to the 
following final expressions 

M-(r, r') = T-^,(r, r') + T-„^(r, r') + T-,,(r, r') 

+5{r - r') J Klliv, v")dv" - 7io(r)9^a,y,,(r) 

-Kl{vy), (13) 

where 

TTf^^A^: r') = do.d'p5{v - r')T-(r), (14) 

= ^ + a;a.) [pi (r, r') - Pi5(r, r')] I r,^^ , (15) 

pi(r,r') is one-particle density-matrix, 

K^liv^v') = no(r)no(r')[.9(r,r') - l\d^,d,Vc{\v - r'|), 

(16) 



(''') .g(r, r') is the pair correlation function defined asi^ 



g(r,r') = ^.^^ /,\ (>5(r - v,)5{v' - r^)), (17) 
no(r)7io(r') 

where ri and Y2 are the coordinates of two distinct elec- 
trons of the system, and T4c(r) is the static xc potential. 



4 



For the scalar f^c, the substitution of Eq. (|13l) into mechanics" approach is to introduce a deformed ground- 
Eq. leads to state wave functionS^ 

/,°?(r,r') = J 3[a;'5;'^(r,r")]T;,^(r")[a::a:'g(r",r')]dr" Vo[u](ri, ...r^,) ^ ^o(ri-u(ri), r^-u(r^)) fj g-V4(r, 

^ z=l 

(20) 

where the factors g ^''^ ensure the correct normalization 
of the deformed state, g^^^ being the Jacobian of the non- 
volume-preserving transformation of coordinates r — > r — 
u(r),i. e. the determinant of the matrix (5^jy—(?^Ujy. Then 
we have 



+ /[a;'g(r,r")]l^S^(r")[a;'g(r",r')]dr" 

[5;:a(r,r")]A--(r",r"')[a-g(r"',r')]dr"dr"', (18) 
where an effective "xc curvature" i^^1{r) is defined as 
f]- (r) = J Klliv, r')dv' - no{v)d,.d.V.4r). (19) 



IV. FORMULATION IN TERMS OF 
CONTINUUM MECHANICS 

In spite of their complicated appearance, the physical 
significance of formulas ([T5|) - ([T^ is quite transparent. In 
fact, the formulas could be derived more quickly start- 
ing from the representation of M^^(r,r') as the second 
functional derivative of the exchange-correlation energy 
functional Exc[v] with respect to virtual displacements, 
Uf_i{r), Uu{r'). The explicit expression for -E'a;c[u] can 
be e asily obtained from the formulas given in Refs. [l^ 
and hJ for the full energy functional E[u] (kinetic plus 
Coulomb interaction plus external potential energy) by 
subtracting the non-interacting kinetic energy functional, 
the Hartree energy functional, and the Kohn-Sham po- 
tential energy. The main idea behind this "continuum 



S^E[u] 



(21) 



u=0 



where £^[u] is the expectation value of the Hamiltonian 
in the deformed ground-state '00 [u]. Similarly, we have 



M^^(r,r') 



S^E^,[u] 



(22) 



u=0 



where, as noted above, i?xc[u] is defined as the difference 
E[u] - Eh [u] - Ts [u] - Vs [u] , where Eh [u] is the Hartree 
energy, Tg [u] is the kinetic energy of the deformed Kohn- 
Sham ground-state (i.e., the non- interacting ground state 
that has the same density as the true ground-state), and 
Vs[u] is the Kohn-Sham potential energy in the deformed 
ground state. 

The above expressions are very helpful in understand- 
ing the physical significance of the terms that appear in 
the evaluation of the double commutator ([S]). To second 
order in u the xc energy functional has the ioimM^ 



E^c[u] = 



+ ^ I drdr'K(r) - u^(r')]A-^(r, r')K(r) - u,{r')] , 

I 



(23) 



where T'^^(r) and K'^l{r, r') arc defined by Eqs. (|15p and 
p6)) respectively. The second term in ((23|) is a typi- 
cal energy of an elastic medium, which is quadratic in 
the deformation gradients (strains).— The correspond- 
ing elastic moduli are determined by the kinetic stress 
tensor T^i/. The third term is an additional nonlocal 
contribution which should be present if the continuum is 
charged and the displacement produces a local polariza- 
tion. Hence the third term corresponds to the xc part 
of the dipolc-dipole interaction energy of an inhomoge- 
neously polarized continuum. Apparently the first two 
terms vanish in the case of a homogeneous (rigid) dis- 
placement since it does not cause deformations and thus 
does not cost any energy. For the rigid displacement 
only the first term in ([23]) survives. This term ensures 



the harmonic potential theorem. It cancels the xc part of 
the KS potential, which guaranties that the rigid motion 
of the whole system is controlled solely by the external 
potential. 

The above expression for E^c can be rewritten in terms 
of the effective curvature il^^{r), defined in Eq. and 
we get 

+ i j dr{n'^^u^u^ + T^l[Au^aUi,a-{d^,Ua){d^Uo)]] , 

(24) 

from which the expression (|13p for M^^ can be straight- 
forwardly recovered by isolating the coefficient of 
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M^(r)u^(r') (this requires some integrations by parts). 

n^\G- G")AI^^„^,„ (q)no ' (G'"- G') , (25) 

V. APPLICATION TO THE PERIODIC CASE 

In the case of a periodic system, the response functions 
and xc kernels become infinite matrices indexed by re- 
ciprocal lattice vectors. The Fourier transform of Eq. (|6]) 

gives where tiq ^(G) is the Fourier transform of l/7io(r) and we 

^ ^ (, imply summation over the repeated reciprocal vectors. 

/gg' Lii^i'iT^) ~ — 9 ^ By Eq. (fT3l) . we can write in the reciprocal space 



MSG',,.(q) = (G + q) • (G' + q)T;,=(G - G') + (G^ + g„)T„-^(G - G')(G. + g.)+ 
(G;, + g^)(G„ + fc)r-(G - G') - [A'aG',^.(q) - K^Go,^,iO)6oG'] 
+(G; - G'piGl - G'MG ~ G")Ke(G" - G'), 
where Tf^^{G) is the Fourier transform of T^l;{r) and K^f^t „^(q) is the double Fourier transform of A'j^^(r,r'). 



(26) 



We can write 



where 



/GG~(q) = 5GG"(q)Afa';-G"'(q)aG"'G'(q), 

A^GG'(q) - (G, + 9^)M-G'.^.(q)(GL + q.) = (g, + q,,)x 

{3(G + q) • (G' + q)T-(G - G') - [A'a'b',^.(q) - K^^^uWSoG'] 
+ {G'^ G'PiG', G'J)no(G - G")V.c(G" - G')} (G', + g.). 

I 



(27) 



VI. LONG-WAVELENGTH LIMIT 

Because of the applications in optics and because of 
the singularities arising in the scalar /^^, the long- wave 
limit (q — 0) requires a separate consideration. 

By the structure of Eq. ([26| , -/1/gg' /it/ (q) have finite 
limit as q — !> 0. By Eq. the same is true for 

/gg' Aii/(q)- We note that in the long- wave limit the 
"head" and "wing" elements of A/qq, reduce to 

M^o,^. = AfoG,.,, = G;:G'>o(G - G")F.e(G"), 

i.e., they are expressed in terms of the ground-state den- 
sity and xc potential only. Equation ()25p . however, shows 
that this is not the case with ./oo.^i^ ^^'^ /G#o,o,/i!/' 
knowledge of which Tqq, and ATgg' iiu needed. 

For the scalar f^'^, it can be, however, noticed that 
Eq. p7p gives a singular long-wave limit. Indeed, since 
5GG'(q) =A/'GG,(q), where 



and the upper row (left column) of A/qg' (q) is zero in 
the q ^ limit, 5GG'(q) is singular in this limit. To 
isolate the singularity, we introduce a matrix 



A/'GG'(q) 



A^GG' (q) 



|G + q||G' + q| 
the inverse of which is regular at q 0. We can write 

•A^GG"(q)Ma';.G"'(q)A^G'"G'(q) 



/gg' (q) 



|G + q||G" + q||G"' + q||G' 



(29) 



Retaining the leading terms in 1 /q only, we have in the 
long- wave limit 



/oT°°(q^0) = 3, 



(30) 



AAGG'(q) = (G + q) • (G' + q)7io(G - G') (28) 



where 
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- ^|-^Go(q)'7MA^oo,^.9.-<"o'(q)+ E [■^GG'(q)A^o"o'(q) + +A^Go(q^^^ 



G'^0 



J2 ^gg'{^)G',M-g'g'',,uG:Mg''o{(i) 



/G#o,Gvo(q ^ 0) 



^ {■A/'Go(q)9MA^oo%.9.A^oG'(q) 



G"5^0 



"-^^-'"G"o,pi/yi' 



E -^GG" (q)G;;MG"G'" ..uGi'^G^'G' (q) 

G",G"'7tO 



r 



A/qq, arc taken at q = and vectors with a hat 
denote unit vectors in the same direction. 

The singularity in q of the "head" and "wings" ele- 
ments of the scalar /^'^ in the long-wave limit has an 
important physical implication: it is directly relevant to 
the description of the excitonic effect in semiconductors 
and insulators by means of TDDFTi^ii^ We point out 
that this singularity is entirely due to the inhomogene- 
ity, since it disappears in the homogeneous case when the 
non-diagonal elements of the matrices are zero (cf. Rcf. 



VII. SUMMARY AND CONCLUSIONS 

In the high frequency limit, we have worked out the dy- 
namical exchange-correlation kernel of an inhomogeneous 
many-body system in terms of the ground-state (or equi- 
librium at finite temperature) properties of the system. 
Both the scalar xc kernel fxc and tensor xc kernel fxa 
relevant to the time-dependent density-functional theory 
and the time-dependent current- density functional the- 
ory, respectively, have been obtained. 

The ground state properties which define the infinite- 
frequency xc kernel have been found to be: (i) Parti- 
cle density; (ii) The kinetic stress tensor; (iii) The pair- 
correlation function; (iv)The exchange-correlation poten- 
tial. The former three can be found, in principle exactly, 
by the method of the quantum Monte-Carlo, while the 
latter can be obtained by the inversion of the Kohn-Sham 
equations after the ground-state density has been found. 



Keeping in view applications to finite and non-periodic 
infinite systems, we have worked out the real-space form 
of the high frequency xc kernels. The real space for- 
mulation is most directly connected with the quantum 
continuum mechanics formulation of Ref. [l3l and clearly 
displays the physical meaning of the anti-adiabatic fxc 
as clastic constants (stiffnesses) of the system. 

In the periodic case and in the long-wave limit, our re- 
sults contain an important singularity in the scalar fxc as 
a function of the wave- vector. This singularity is known 
to be a manifestation of the ultra non-locality of fxc in 
space and its presence is crucial for reproducing the ex- 
citonic effect in semiconductors and insulators by means 
of TDDFT. 

We expect the results of this work to facilitate the con- 
struction of exchange-correlation kernels of TDDFT and 
TDCDFT, accurate enough to account for the many- 
body effects in the linear response theory with appli- 
cations to optics, electron energy-loss spectroscopy, and 
other fields. In particular, approximate formulas may be 
constructed for the frequency dependence of the xc ker- 
nels, interpolating between the adiabatic and the anti- 
adiabatic limits. In any case, it remains an interesting 
question to fully understand the relation between the adi- 
abatic and the anti-adiabatic approximation to the xc 
kernels. They are formally very similar, instantaneous 
functionals of displacement and density respectively, but 
the physics they embody is quite different since the sys- 
tem is assumed to be fully relaxed to the instantaneous 
ground-state in the former case, while, in the latter, the 
elastic deformation dictates the instantaneous state of 
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the system. 
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Appendix A: Derivation of M^^;(r, r') 

In this Appendix, we outline the derivation of the ex- 
pectation values of the commutators relevant to the eval- 
uation of the moment by Eq. ([S]) . Since the system is as- 
sumed paramagnetic, the equilibrium expectation value 
of the current is set to zero throughout. 

(i) Kinetic energy part. A straightforward evaluation 
leads to 



Introducing a tensor 



we express it as 

where pi{r, r') is the one particle density matrix. For the 
kinetic energy part we, therefore, have 

m^([j,ir), [T,Mr'm = r')T^.(r)+ 



+ ^a;a,a;9„(5(r-r')no(r). 



(Al) 



(ii) The external potential part is easily evaluated to 



no{r)d^d{r-r')d,V„{r). (A2) 
(iii) For the interaction part, we need the two types of 
commutators 



m2 [n^(r), [Vc{\ri - r^l), ji,{r')]] = 6{r - n)d^S{r - r')dlVc{\r' - r^l), 
b'l^r), [Uc(|ri - ral), j2.(r')]] = '^(r - n)6(r' - r2)d^dlVc{\r ~ r'|). 

Then the interaction part can be rewritten in terms of the pair correlation function g(r, r') of Eq. (jl7p as 



i5^yc(|r„-r„|),j.(r') 



n^ra 



= no{r)x 



[g(r,r") - l]d^S{r - r')9,>o(r")l^c(|r' - r"\)dr" + no(r')[.9(r, r') - l]d^dlVc{\r - r'|) 

-d^dir - r')d'Mr') + noir')d^dlVci\r - r'|)] , 

I 



(A3) 



where 



ton's law 



Vnir) - / noir')Vc{\r^r')\)dr' 



is Hartree potential. 

(iv) A considerable simplification comes from the New- 



9,r,^(r)+no(r)a^yo(r) 

+no(r) J 5(r, r')d/^^^dr" = 0, (A4) 

the latter immediately following from the obvious rela- 
tion {[j^ir),Ho]) = 0. 
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(v) 
1 



Mp^(r,r') is obtained by combining Eqs. (|Aip . 
and Then M^^(r,r') is constructed 

from A/p^(r, r') by replacing Vb with Vs ^Vo + Vh + Vxc 
and setting the interaction term to zero. We finally arrive 
at Eq. ([TS]) by the use of Eq. 



Appendix B: The anti-adiabatic limit of the xc 
kernel for the homogeneous electron gas 

The tensor xc kernel of a homogeneous electron gas has 
the form^ 

(Bl) 

where q is the unit vector in the direction of q. 

The longitudinal (L) component of the kernel coincides 
with the scalar xc kernel. In the high-frequency limit, 
by using the f-sum rule and the third-moment sum rule 
for the density-density response functions (or the first- 
moment sum rule for the current-current response func- 
tions) for the interacting and non-interacting gas,— one 
readily obtains 



fxcLil^Oo) 



(B2) 



where 

and 
with 



:ar: / dk [S{k,r,) - 1] f (k/q) (B3) 



(B4) 




FIG. 1: High-frequency limit of the longitudinal xc-kernel for 
the homogenous electron gas plotted as function of q for differ- 
ent values of rs. Asymptotically the only finite contribution 
is due to the term containing the correlation kinetic energy. 




X 

Y 



2.T 



■In 



1 + 1/x 



1 - 1/x 



tc{rs) is the correlation kinetic energy, and S{k, r^) is the 
static structure factor. We emphasize that the contribu- 

(2) 

tion due to the correlation kinetic energy, f^jJi^iq^oo), is 
independent of q. In the expressions above the atomic 
units in which h = = m = 1 are used and the wave 
vectors are expressed in units of the Fermi wave vector 
qp = a/rs, with a = (97r/4)^/^. 

In the same limit, the transverse component is given 
by 



fxcrili 



(B6) 



For the calculation of S{k,rs) we have employed the 
analytic static structure factor provided in Ref. [l^ and 
the correlation kinetic energy has been calculated from 
the parametrized correlation energy of Appendix B of 
the same reference, by making use of the virial theo- 
rem. The resulting quantities f'^chi^^ fxcxil^ °o) and 



FIG. 2: Here the quantity in Fig. ([T]) is plotted by neglecting 
/fi(g,oo) [Eq. mi 



/ic^((7, oo) arc plotted in Fig. Fig. © and Fig. © 
respectively. Among other facts, these plots show that 
the contribution due to the correlation kinetic energy is 
very important at any g - a fact already pointed out 
by Iwamoto and Grosa^^ for the longitudinal component 

at g = 0. Moreover, it is clear that fxcL^I^'^) '-^^ 
be identified as the q ^ oo limit of f^^^{q, oo). Over- 
all, /^cL(T)(9' °'^) bas a strong dependence on both q 
and Ts- From Fig. ([T]) and Fig. ([3]), it is evident that 
/kcl(9'Oo) ^^'^ fxcri-l^'^) have very different ranges of 
values. We observe that while f!^^rp{q,oo) is always pos- 
itive, /^^^(g, oo) changes sign at about a fixed value of 
q/kp for all rg'- this suggests that the r^-dependence of 
/^'^^(g, oo), with q expressed in units of kp may be an 



overall scale factor. Clearly, the fact that /' 



cL(T) 



(q, oo) 



go to a constant in the large q limit implies that the 
corresponding local field factors diverge as g^. This is 



9 




also observed at finite frequency, and, indeed all the 
way to zero frequency, as first noticed by HolasSS and 
then confirmed by quantum Monte Carlo^ and diagram- 
matic calculations 1^ In particular, the exact relation 
oo, oo) = —'ifxchil ~^ is well satisfied. 



FIG. 3: High-frequency limit of the transverse xc-kernel for 
the homogenous electron gas plotted as function of q for differ- 
ent values of rs- Asymptotically the only finite contribution 
is due to the term containing the correlation kinetic energy. 
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